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Spin- 1/2 two-legged ladders respecting inter-leg exchange symmetry o and spin rotation symmetry 
D2 have new symmetry protected topological (SPT) phases which are different from the Haldane 
phase. Three of the new SPT phases are tx,ty,tz, which all have symmetry protected two-fold 
degenerate edge states on each end of an open chain. However, the edge states in different phases 
have different response to magnetic field. For example, the edge states in the tz phase will be split 
by the magnetic field along the z-direction, but not by the fields in the x- and y-directions. We 
give the Hamiltonian that realizes each SPT phase and demonstrate a proof-of-principle quantum 
simulation scheme for Hamiltonians of the to and tz phases based on the coupled-QED-cavity ladder. 
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I. INTRODUCTION 

Symmetry protected topological (SPT) phases are 
formed by gapped short-range-entangled quantum states 
that do not break any symmetry.^ Contrary to the trivial 
case, quantum states in the non-trivial SPT phases can- 
not be transformed into direct product states via local 
unitary transformations which commute with the sym- 
metry group. Meanwhile, two different states belong to 
the same SPT phase if and only if they can be trans- 
formed into each other by symmetric local unitary trans- 
formations.^ A nontrivial SPT phase is different from the 
trivial SPT phase because of the existence of non-trivial 
edge states on open boundaries. This non-trivial prop- 
erty is protected by symmetry, because once the symme- 
try is removed, the SPT phases can be smoothly con- 
nected to the trivial phase without phase transitions.'^ 
The well known Haldane phase"* is an example of SPT 
phase in 1-dimension (ID). Topological insulators^"*" are 
examples of SPT phases in higher dimensions. 

The bosonic SPT phases are classified by projective 
representations (which describe the edge states at open 
boundaries or at the positions of impurities**"*'') of the 
symmetry group in ID,*^ and by group cohomology the- 
ory in higher dimensions.*^ We also have a systematic 
understanding of free fermion SPT phases*® and some in- 
teracting fermion SPT phases. *^"^° Using those general 
results, fourteen new ID SPT phases protected by 
spin rotation and time reversal symmetry are proposed 
in Ref. 21. 

In this work, we will discuss two-legged spin-1/2 lad- 
der models with two-body anisotropic Heisenberg inter- 
actions respecting D2 x a symmetry. Here a is the inter- 
chain exchanging symmetry and D2 — {E, R^, Ry, Rz}, 
where E is the identity, R^ {Ry, Rz) is a 180° rotation of 
the spin along x {y, z) direction. The symmetry D2 x a 
protects seven non-trivial SPT phases. Four of them, tg, 
tx, ty, tz, can be realized in spin-1/2 ladder models. The 
to phase is the Haldane phase, and the tx, ty, tz phases 



are new because of their different edge states. We pro- 
vide a simple two-body Hamiltonian for each SPT phase, 
and study the phase transitions between these phases. 
We also discuss possible physical realizations of the SPT 
Hamiltonians and demonstrate a proof-of-principle im- 
plementation scheme based on coupled quantum electro- 
dynamics (QED) cavity ladder. 

The paper is organized as follows: in Sec. H, we in- 
troduce the projective representations of the underlying 
symmetry group of the two-legged spin-1/2 ladder and 
discuss the possible SPT phases. We then numerically 
study the phase diagram and the phase transitions of the 
spin-ladder Hamiltonians in Sec. HL In Sec. IV, we dis- 
cuss the physical realizations of the SPT Hamiltonians 
and demonstrate a proof-of-principle quantum simula- 
tion scheme based on QED cavity ladder. Finally, we 
summarize in Sec. V. 



II. PROJECTIVE REPRESENTATIONS OF THE 
SYMMETRY GROUP AND SPT PHASES 

The symmetry group D2 x a for our two-legged spin- 
1/2 ladder is Abelian. All its eight representations are 
1-dimensional. The following two-spin states on a rung 
form four different ID representations of the symme- 
try group: |0,0) = ^(| ti;2) - I ;it2)), = 

Ui2) - 1 tit2», |i,2;> = ;i;2> + 1 tit2)), 

and |l,z) = ^(1 tii2) + I iit2))- The subscripts 1,2 
label the different spins on the same rung. The group 
D2 X cr has eight projective representations (see Tab. I), 
which describe eight different SPT phases^^ of spin lad- 
der models.*"* One of the projective representations is 
one-dimensional and trivial. The other seven non-trivial 
ones are two-dimensional, which describe the seven kinds 
of two-fold degenerate edge states of the seven non-trivial 
SPT phases. 

The degeneracy at the edge in each non-trivial SPT 
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TABLE I: Projective representations and the corresponding 
SPT phases for D2 x a group. The active operators can spht 
the degeneracy of the ground states. The operator 0± = 
Oi ±02, where 1,2 are the labels of the two spins at a rung, 
and SS- means Si,i ■ Sia+i — 824 ■ S2,i+i. 
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"If the system has translational symmetry, there are four different 
rung-singlet phases: rung-|0, 0) phase, rung-|l,x) phase, rung-|l,iy) 
phase and rung-|l,2). If the system does not have translational 
symmetry, then there is no difference between the four rung-singlet 
phases and there will be only one rung-singlet phase. 



phase is protected by the symmetry. To lift such a degen- 
eracy, we need to add perturbations to break the D2 x a 
symmetry. The operators that split the edge degeneracy 
are called active operators. The active operators are 
different in different SPT phases, which provide experi- 
mental methods to probe different SPT phases. 

We will mainly study five of the eight SPT phases 
which can be realized in two-legged spin- 1/2 ladders: one 
trivial SPT phase (corresponding to £"0) and four non- 
trivial ones to, tx, ty, tz (corresponding to E^,E^,Eq, Ej 
respectively). The to phase is equivalent to the spin-1 
Haldane phase, where the edge states behave like a spin- 
1/2 spin and will be polarized by a magnetic field in arbi- 
trary direction. However, in the t^ phase, the edge states 
will only respond to the z component of a homogeneous 
magnetic field ('homogeneous' means that the field takes 
the same value at the two sites of a rung) . That is to say, 
the ground state degeneracy will be lifted by but not 
by Bj. or By (see Appendix A) . The t^, and ty phases are 
defined similarly. From the response of the edge states 
to external magnetic field, we can distinguish the four 
phases io, ^x, iy, ^z-^^ 

The other three SPT phases can be realized by stacking 
two of the above ones: the SPT phases corresponding to 
El, E2, E4 can be realized by stacking t^ and ty, ty and 
tz, tx and tz, respectively. The to phase can also be 
realized by stacking three phases tx,ty,tz. 

The nontrivial SPT phases tQ,tx,ty,tz can be realized 
in two-legged spin-1/2 ladders. The Hamiltonian that 
realize these phases is simply given a.s H — Hl + Ht, 
where Hl is the interaction along the leg (or longitudinal 
direction) and Hx is the interaction along the rung (or 
transverse direction). We assume takes the following 



form, 

— [('^'5'i,i'S'i_i+i + S\^Sf ^^i + <S'f i5'J^i+i) 

i 

+ {^^2,1^2,1+1 + ^2A^2,i+l + ^2A^2,i+l)] J 

where J > 0, S ^ I, and 5™j (m = x,y,z) is the spin 
operator at the jth leg and ■ith rung. 

By tuning the interaction Ht, we can obtain four SPT 
phases tQ,tx,ty,tz. For example, the following model 

Ho = i/i + A^Si,, -82,,;. (1) 

i 

can realize the to phase. When A > 0, the rung-singlet 
state (|0, 0)) is lower in energy and the system falls into 
the rung-|0, 0) phase, which is a trivial SPT phase. When 
A < 0, the rung-triplet x) ,\l,y) z)) are lower in 
energy, and effectively we obtain a spin-1 anisotropic 
Heisenberg model, which belongs to the Haldane phase 
to- 

By partially flipping the sign of interactions along the 
rung, we obtain the Hamiltonian for the tz phase: 

Hz = Hl + X Y.{-Sl,Sl, - S^Sl^ + Sl,Sl,).{2) 

i 

When A > 0, it falls into a trivial SPT phase which corre- 
sponds to the rung \l,z) product state. When A < 0, the 
low energy degrees of freedom in each rung are given by 
three states \l,x),\l,y),\Q,Q), and the resultant model 
belongs to the tz phase. Note that the Hamiltonians 
(1) and (2) (with the respective ground states to and 
tz) can be transformed into each other by the unitary 
transformation: Ui{tt) — exp (ztt S'f j) . However, 
this does not mean that to and tz belong to the same 
phase since Ui{t:) does not commute with the symmetry 
group D2 X a. Applying similar arguments, we may have 
Hx = Hl + \Y.^{Sl,Sl, - Sl,Sl, - Sl^Sl,) for the tx 
phase, and Hy - i^L+A 

for the ty phase. Tab. I shows that different SPT phases 
have different active operators. For example, in the to 
phase, the edge degeneracy can be lifted by either 
or S\ or S^. This means that the edge states can be 
polarized by a homogeneous magnetic field along any di- 
rection. In the tz phase, S'^ and S\_ are not active opera- 
tors, indicating that a weak homogeneous magnetic field 
in the x-y plane will not split the edge degeneracy. These 
properties can be verified by a finite-size exact diagonal- 
ization study of the Hamiltonian (1) and (2).^^ 

III. PHASE TRANSITIONS AND THE PHASE 
DIAGRAM 

To study the phase transitions, we consider the follow- 
ing model which contains both to and tz phases, 

h = Hl + Y.\^.y{si,si, + si^si;) + XzSl^Sl,]. (3) 
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FIG. 1: (Color online) The phase diagram of the Xz-Xxy plane 
with S = 0.9. The point (\xy,\z) = (0,0) is a multi-critical 
point. Near this point, the phase boundaries are hard to 
determine. The dashed line means that the phase bound- 
aries are not accurate. The two 'XY' phases are gapped and 
magnetically ordered in y-direction because 5 7^ 1 introduces 
anisotropy in the x-y plane. The fo/iz SPT phase locates be- 
tween an Ising-type ordered phase and an XY-type ordered 
phase. 
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FIG. 2: (Color online) Phase transitions with 5 = 0.9, Az = 
—0.5 J. The order parameters are defined as M™ = \l^ST ± 
S™)]. The green line with asterisks shows the information of 
entanglement spectrum Ap = pi ^ p2, where pi and p2 are 
the two maximum Schmidt eigenvalues of the ladder, Ap = 
means that the entanglement spectrum is doubly degenerate. 
The SPT phases are characterized by doubly degeneracy in 
the entanglement spectrum and vanishing of all magnetic or- 
ders. 



Since the Xxy-^z plane phase diagrams of above model 
are quite different for 5^1 and 5 = 1, we will discuss 
them separately. 



A. the case 5 7^ 1 

1. The phase daigram 

We only consider the case that anisotropy is weak, so 
we set S = 0.9 for simplicity. The phase diagram Fig. 1 
(with 5 = 0.9) is obtained via infinite time-evolving block 
decimation (iTEBD) algorithm. The phase diagram is 
symmetric along the line Xxy = 0, because the model with 
—Xxy can be obtained from the one with X^y by the uni- 
tary transformation Uii^n). The origin {Xxy,Xz) — (0,0) 
is a multi-critical point linking all the phases. On the 
upper half plane Xz > 0, there are only three phases. 
The lower half plane is more interesting. The limit 
Xz —00 corresponds to the Ising Stripe-Neel phase, 
while Xxy — ±00 corresponds to the rung-|0, 0)/rung- 
|1, z) phase. 

In the intermediate region, we have two XY-like phases 
and two SPT phases: XY-stripe Neel phase and XY-Neel 
phase. For (5 7^ 1, the exchange interactions along the 
legs are anisotropic in the x-y plane, consequently the 
two XY-like phases have a finite excitation gap and are 
ordered in y-direction if (5 < 1 (or ordered in x-direction if 
S > I). The to/tz phase is located between the XY-stripe 
Neel/XY-Neel phase and the Ising stripe Neel phase. The 
phase diagram at Xz = — 0.5J is also shown in Fig. 2, 
where we illustrate the symmetry breaking orders and the 
entanglement spectrum (Ap = pi— p2 is the difference be- 



tween the two biggest Schmidt eigenvalues) in each phase. 
For 5 = 1, the two XY phases vanish, and the SPT phase 
to/tz touches the trivial phase rung-|l, z)/rung-|0, 0) di- 
rectly (we will study this case in detail later). 



2. Semiclassical explanation of the phase diagram 

Notice that both of the SPT phases are sandwiched 
by two ordered phases. This suggests that they originate 
from quantum fluctuations caused by the competition be- 
tween the different classical orders. To understand this 
point better, it will be interesting to compare the phase 
diagram Fig. 1 with that from the semiclassical approach 
in which the ground state is approximated by a direct 
product state. A semiclassical phase diagram (Fig. 3 
with 5 = 0.9, Xz = — 0.5J) is obtained by minimizing 
the energy of the following trial wave function, 

|V')sc = ]^(ai'/'l,2! + a202,2i + a303.2i + 0404, 2i) 

i 

<8'(a201,2i+l + ai02,2j+l - a3(t)3^2i+l + a4'?!'4,2i+l), 

where 0i = -^{\l,x) +i\l,y)), 02 = -^i\l,x)-i\l,y)), 
03 = |l,z), 04 = |0, 0) are four bases of each rung, and 
01,02,03,04 are four trial parameters. 

Comparing the semiclassical phase diagram with 
Fig. 2, we find that the two phase diagrams are simi- 
lar except for the absence of two SPT phases in Fig. 3. 
Importantly, the locations of the quantum SPT phases 
are close to the phase boundaries between the Ising or- 
dered phase {M^ ^ 0) and the xy-planar ordered phases 
(M^ 7^ 0) of the semiclassical phase diagram. Similar 



situations occur in spin-1 XXZ chain, where a Haldane 
phase exists between two ordered phases. ^'^^ This sug- 
gests that we can roughly obtain a quantum phase dia- 
gram with a semiclassical approach by 'inserting' a SPT 
phase at the phase boundary of two different classically 
ordered phases. 

Furthermore, the semiclassical picture even indicates 
some important information of the quantum SPT phases. 
For instance, it tells us why the tp phase is different from 
the tz phase. The classical phase boundary between the 
two ordered phases is located at \X%y\ = |Az|. At the 
point \%y = Az < 0, the states have the same lowest en- 
ergy if the spins are antiferromagnetically ordered along 
the leg and are parallel along the rung polarizing in y-z 
plane (namely, the classical ground states are highly de- 
generate). Quantum fluctuations (which is enhanced by 
the degeneracy) will drive the system into the higher en- 
ergy states (e.g. some spins are pointing off the y-z plane) 
with a certain weight. As a consequence, the ground state 
is short-range correlated with a finite excitation gap. Fur- 
thermore, edge states exist at each boundary. The edge 
staes can be considered as an effective spin whose mag- 
netic momentum is half of the total momentum of the two 
spins at a rung. Since the two spins in the same rung are 
parallel, the edge states carry free magnetic moment and 
respond to magnetic field along arbitrary directions. On 
the other hand, at the point X^y — — Az > 0, the lowest 
energy classical states are those that the two spins are 
parallel along the z direction and anti-parallel in the x-y 
plane (see Fig. 4b). Quantum fluctuations around the 
sphere drive this state to the phase. In the phase, 
the effective edge spins have no net magnetic moment in 
the x-y plane, so they will not respond to the magnetic 
field in the x-y plane. Similar things happen in the tx 
and ty phases. Above arguments are applicable to other 
systems with different symmetry groups, and provide a 
guidance to seek SPT phases and analyze their physical 
properties. 
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FIG. 3: (Color online) The order parameters and in 
the semiclassical approach. Compared to Fig. 2, the two SPT 
phases are missing. 
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FIG. 4: (Color online) (a) At \%y = A^ < 0, the classical 
states on the sphere are close in energy. Each direction on 
the sphere defines a classical state in which two spins at the 
same rung are parallel and pointing in that direction, (b) At 
Xcy = —X^ > 0, the classical states on the sphere are also 
nearly degenerate. But now the two spins at the same rung 
are antiparallel in the x-y plane and parallel ofi' the plane. 



B. 



the case 5=1 



1. Finite size effect in numerical method 

When establishing the phase diagram Fig. 1 and Fig. 2, 
we have approximated the ground state by a matrix 
product state (MPS), which is obtained through iTEBD 
method.^'* The matrices in the MPS have a finite di- 
mension _D, which introduces a cutoff to the number of 
Schmidt eigenvalues in the entanglement Hilbert space. 
While the energy of the actual ground state can be es- 
timated using a finite-D MPS with high accuracy, the 
order parameters are usually overestimated due to the fi- 
nite dimension of the matrix. Furthermore, in some cases 
the order parameters (which are finite if D is finite) even 
vanish in the infinite D limit. So a scaling of the order 
parameters with respect to the dimension D is necessary 
in order to obtain the correct phase diagram. In the fol- 
lowing we will illustrate, via finite size scaling, that the 
two XY phases disappear \i 5 — 1. 

The phase diagram in Fig. 5(a) with 5 — l^Xz — —0.5 J 
is obtained by setting D — 30. To evaluate the order 
parameters of each phase in the infinite D limit, we per- 
form numerical calculations with different D's, as shown 
in Fig. 6. Performing finite size scaling and extrapo- 
lating to infinite -D, we find that the maximum value 
of M^'' in the XY-strip Neel phase (and the same for 
Mj' in the XY-Ncel phase) vanishes as a power law of 
£)-o. 24142^ In contrast, the maximum value of in the 
Ising stripe Neel phase approaches a finite number ex- 
ponentially fast. This shows that when 5 = 1 the two 
XY phases in Fig. 5(a) are due to finite size effect (notice 
that decays very slowly with increasing D, so any 
finite D will give incorrect results) and will disappear in 
the phase diagram in the limit I? — > oo (see Fig. 5(b)). 

To illustrate that the XY phases will not disappear if 
(5 ^ 1, we perform a scaling for the order parameter 
with £) for 5 = 0.9. It turns out that still varies as 
a power law of £)-3-803^ when shifted by a constant Mg 
(See Fig. 7). The fact that Mq ^ Q shows that Ml in the 



5 



XY-strip Neel phase (and the same for in the XY 
Neel phase) is finite in the hmit Z? — > cx), and the phase 
diagram at finite D are quahtatively correct if (5 7^ 1. 



2. Vanishing of the XY phases and (7(1) symmetry 
protected phase transitions 

PhysicaUy, the disappearance of the two XY phases at 
i5 = 1 can be explained by quantum fluctuations. 

When 6 = 1, the Hamihonian (3) has an enhanced 
symmetry [C/(l) x Z2] x a, where the C/(l) subgroup 
means rotation of the spins along z-axis, Z2 is gener- 
ated by a rotation of the spin for n around a;-axis. Since 
the continuous symmetry U{1) will never spontaneously 
break in ID, if the XY phases exist, they will be quasi- 
long ranged ordered in x-y plane and gapless. 

However, strong quantum fluctuations gap out these 
states and drive them into the SPT phases. From the 
semiclassical approach, the enhanced symmetry results 
in a larger degeneracy of the semiclassical ground states 
near X^y ~ . The enlarged degeneracy of the semiclas- 
sical ground states enhances the quantum fluctuations 
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FIG. 5: (Color online) (a)Phase diagram obtained with D = 
30 for 5=1, Xz = —0.5 J. The order parameters are defined 
as = \{St + 51)1 and M^" = ^ {Sf ± S^)^ + {Sf ± S'^y ; 
(b) Phase diagram of 5 = 1 in the limit D — >■ cxa. 



«5 




0.611 



0.610 



0.609 



0.608 



0.04 0.06 
0.10798D 



0.08 



e 

(b) 



FIG. 6: (Color online) (a)Power law scaling of the maximum 
value of the order parameter M^" with 6 = 1. The intercept 
shows that M^^ is vanishing in power law with the dimension 
D. Insert I shows the data of Af^" vs. dimension D. Insert 
II is a log-log fit of the data; (b) Exponential scaling of the 
maximum value of the order parameter with 5=1. The 
intercept shows that A/+ is finite at infinite D. Insert I shows 
the data of vs. dimension D. Insert II is a log fit of the 
data. 



and gap out all the states (except the states at the criti- 
cal points). Consequently the XY phases disappear, and 
the nontrivial SPT phase to/tz touches the trivial phase 
rung-|l, z)/rung-|0, 0) directly. That is to say, the di- 
rect transition between to/tz and rung-|l, z)/rung-|0, 0) 
is protected by the continuous U{1) symmetry. 

Notice that similar situations also happen in 5* = 1 
chains,^ where the Haldane phase and the trivial phase 
(in analogy to the rung-phases) are separated by a Z2 
symmetry breaking phase when the system does not have 
a continuous symmetry [such as U{1) spin rotational 
symmetry] . 
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FIG. 7: (Color online) Power law scaling of the maximum 
value of the order parameter Af^ with S — 0.9. The intercept 
shows that is finite at infinite D. Insert (I) shows the 
data of VS. dimension D. Insert (II) is a log-log fitting 
of the data. 



IV. PHYSICAL REALIZATION AND 
QUANTUM SIMULATION 

Here we propose possible realizations of the SPT 
phases in two-legged ladder Mott systems.^® The Hamil- 
tonian (3) can be considered as two Heisenberg chains 
coupled with magnetic dipole-dipole interaction and ex- 
change interaction Ht = J2i[v{^i,i ' 82,1) + 7(Si,i • ~ 
35';f jS*! j)] with Az = '7 — 27 and A^^j, — rj + j. However, in 
real materials, the magnetic dipole-dipole interaction is 
too weak to support the SPT phases. On the other hand, 
axial anisotropy interaction Dz{Sf + 82)^ = 2Dz{SfS2 + 
1) can yield the Hamiltonian (3) by Ht = J2i[v{^ii ' 
S2,^) + D,{Sl, + Sl,)^] with A, = 2D, + 7j, X^y =\. 
Note that a weak effective spin-1 axial anisotropy term 
exists in the material (C5Hi2N)2CuBr4.^^ If this term is 
negative and is strong enough in some material, then the 
tz phase will be realized. 

On the other hand, due to the nearest-neighbor-only 
interactions, it is tempting to consider the possibility 
of simulating these Hamiltonians in a non-condensed- 
matter setting. Here we provide a proof-of-principlc 
implementation scheme for the Hamiltonians for the 
phase (1) and for the phase (2) based on coupled- 
harmonic-oscillator array. Here we only consider the 
isotropic case with i5 = 1. Note that it is also possible to 
implement the more general anisotropic cases with addi- 
tional Raman lasers along the ladder. Our scheme may 
be extended to systems including solid-spins interacting 
with arrays of coupled transmission line resonators, or 
ions in a Coulomb crystal. As a concrete example, we 
illustrate the scheme using a coupled-QED-cavity ladder 
(see Fig. 8(a)), where the quantized cavity fields cou- 
ple to their nearest neighbors in the longitudinal (L) and 
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FIG. 8: (Color online) (a) Illustration of the coupled-cavity 
ladder. A single atom is held at the intersection of two cav- 
ity modes that are coupled across both the longitudinal (L) 
and the transverse (T) directions, with coupling rates and 
, respectively, (b) Schematic of the coupling scheme for 
the atom inside the cavity. The external lasers are applied 
perpendicularly to the plane of the ladder. The subscripts 
J, k of the parameters indicate the atom is on the jth site 
transversally and fcth site longitudinally. 



the transverse (T) directions via photon hopping.^® The 
spin degrees of freedom on each site are encoded in the 
internal states of a single atom trapped inside the cavity. 

In Fig. 8(b), we show the coupling scheme for a mini- 
mal model with three internal states, with (| tj.fc)? I 4-j,fc)) 
corresponding to hyperfine states in the ground state 
manifold and \ej^k) corresponding to an electronically ex- 
cited state. Two far-detuned external laser fields with 
Rabi frequencies f^^^ = O'^e'^^'^ {^l = L, T) couple |tj,fc) 
and \ej,k), while two far-detuned cavity modes with Rabi 
frequencies couple \ej^k) and |4.j,fe), with all the de- 



tunings satisfying A. 



|Al^-A'f|,|A5'| 



lA^I 



C*, ri^, IA2 — A^* !. Hence, we have two independent Ra- 
man paths in the longitudinal and in the transverse di- 
rections, and a virtually populated excited state \ej.k)- 
Finally, the hyperfine states are coupled using a res- 
onant radio-frequency (r.f.) dressing field with Rabi 
frequency JlJ^j. = Jl'^^e"^^ '' . One can show that all 
the cavity modes are virtually excited under the con- 
dition lA^-A^'l ~ » g>'/V2N,vt', where 5^ = 
Gt^nt'iA^^ + A^)/2A5'A^, is the coupling rate be- 
tween cavities along the direction /x, and N is the num- 
ber of sites in the longitudinal direction. After adi- 
abatic elimination of the excited states and the cav- 
ity modes, the resulting effective Hamiltonian gives rise 
to effective spin-spin interactions between the neighbor- 
ing sites, with the pseudo-spins in the effective Hamil- 
tonian given by the superpositions of the hyperfine 
states: | t;-,^) = {e^^-"] tjj.) + \ hk)) /V2, \ i'^J = 
(-|tM-)+e-'^-Hkfc))/V2. 

Depending on the magnitude and the relative phases 
of the coupling fields, the effective Hamiltonian can sup- 
port either the to phase or the phase. In particular. 



with fl^ 



0.75|A^ 



A^^l and 



^j,k+i 



qL 



the Hamiltonian for the Iq phase when 0jf, — (pj^k 



TT, we get 



= 
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for arbitrary j, fc; and we get the Hamiltonian for the 



tz phase when 



h,k 
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TT and 



ipj^k = Vj,k+i for arbitrary j, k. These conditions can 
be achieved by modulating the relative phases of the ex- 
ternal coupling fields so that they are either or tt be- 
tween neighboring sites. The interaction rates J and A in 
Hamiltonian (1,2) are given as 

A — —lA4:{g^^^v^/(p.'^^Y^. For typical experimental 
parameters'^O; - lOOMHz, C - lOOMHz, |Af| ~ 
IGHz {i = 1,2), - lOMHz, VC^ - lOOMHz, we have 
J ^ O.llMHz, with the magnitude of A/J widely tunable 
by adjusting the ratio between Q,^ and $7"^. 

The single-site addressability of the atoms in the 
cavity ladder allows much freedom in probing the 
properties of the system. For the SPT phases in which 
we are interested, we need to measure the response of 
the edge states to the magnetic field. As the pseudo- 
spins are related to the hyperfine states via a rotation, 
we may implement an effective magnetic field on the 
pseudo-spins by imposing the appropriately rotated 
operations on the edge states. For instance, to apply an 
effective magnetic field in the z-direction, we need to 
apply a rotation on the hyperfine states (assuming 
V'i.fc = 0), which can be achieved for example by adding 
a resonant r.f. field between the hyperfine states, 
with the effective Rabi frequency corresponding to the 
magnitude of the effective magnetic field. To probe the 
response of the edge states to the effective field, we 
may measure the polarization of the ground state of the 
system. Alternatively, as the effective magnetic field 
may induce an energy splitting between the spin states, 
we can measure the existence of energy splitting as a 
response of the SPT phases to the external field. As an 
example, for the phase under an effective magnetic 
field along z, we need to implement the following steps 
on the two edge sites at one end of the open boundaries: 
(i) adiabatically turn on resonant r.f. fields between 
the hyperfine states so that the degeneracy of the edge 
states is lifted and an energy splitting appears between 
I todgo) and I ledge); (n) in the presence of the effective 
magnetic field, apply an effective resonant coupling 
fields between the pseudo-spins, which corresponds to 
two-photon detuned Raman fields between the hyperfine 
states with the Stark shifts equal to the Rabi-frequency 
of the effective resonant coupling fields between the 
pseudo-spins; (iii) after some time of evolution, rapidly 
turn off all the coupling fields, then apply a 7r/2-pulse 
on the hyperfine states, so that pseudo-spin population 
is projected onto that of the hyperfine states; (iv) 
the population of the hyperfine state can be probed 
for example by high-fidelity hyperfine state readout 
technique based on cavity-enhanced fluorescence."^^ If 
the edge states are responsive to the effective magnetic 
field, one will observe Rabi-oscillations in the measured 
fluorescence. Measuring the magnetic field response in 
all three spatial directions will allow us to establish the 
signature of the tz phase. 
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In summary, we have studied nontrivial quantum 
phases to (the Haldane phase) and tz(a new SPT phase 
different from the Haldane phase) protected by D2 x ct 
symmetry in a spin- 1/2 ladder model. The model has 
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We then provided a semiclassical understanding of the 
physical properties of the SPT phases, and discussed the 
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systems. 
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Appendix A: Constructing the Hamiltonians with 
SPT phases 

In this Appendix, we demonstrate the method^^ by 
which we obtain the Hamiltonian Hq, H^, Hy and in 
the main text. The procedure contains three steps: (i) 
construct a matrix product state (MPS) wave function 
with given edge states which are described by the pro- 
jective representations; (ii) construct the parent Hamil- 
tonian for the MPS using projection operators; (iii) sim- 
plify the parent Hamiltonian by adiabatic deformations. 
In the first step, we need to know the projective and lin- 
ear representations of the symmetry group, together with 
the Clebsch-Gordan (CG) coefficients. In the following, 
we will first present necessary information and then dis- 
cuss the construction of the Hamiltonian step-by-step. 

The four projective representations of D2 x a corre- 
sponding to the to,tx:,ty,tz phases are given in Tab. II 
(we only provide the representation matrices for the gen- 
erators) , and the eight linear representations are listed in 
Tab. III. 



TABLE II: Four nontrivial SPT phases in 5 = 1/2 spin lad- 
ders respecting D2 x a symmetry. 





M{Rz) 


M{Rx) 


M{a) 


active operators 


E'i (to) 


ioz 




I 






ioz 


Ox 


iox 






ioz 




Oz 


(5^,5^,5^) 


E? (ty) 


ioz 


iax 


irjy 


(5^,5^,5^) 
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TABLE III: Linear representations of D2 x f. The four bases are defined as |0,0) = ^(| tiiz) — | iit2)), 11,2^) ~ ^(| iii2 
) - I tit2)), = ^(1 iii2) + I tit2)), and 11,2;) = ^(1 tii2) + I iit2)), where the subscripts 1,2 label the different spins 

on the same rung. The operator 5+ = 5™ + 5*™ (m = x, y, z) is even under inter-chain reflection, and 5™ = ST — S*™ is odd 
under the reflection. 





R. 


Rx 




bases 


operators 


Ag 


1 


1 


1 






Big 
B2g 
Big 


1 

-1 
-1 


-1 
-1 
1 


1 
1 
1 


|i,y> 

\l,x) 


s+ 

SI 
SI 


Au 


1 


1 


-1 


|0,0) 




Blu 


1 


-1 


-1 




SI 


B2u 


-1 


-1 


-1 




SI 


Bsu 


-1 


1 


-1 




5^ 





<)E3 








3 Bsg, 


=cr^, C^is 






C^^'s /; 




^E5 








3 Bsu, 




= /, C^is = 


i<Jy, 






^Eq 








B Blu, 






^I, 


C"'^ = a,; 


Er6, 


^Er 



















The Hilbert space of the direct product of two projective representations can be reduced to a direct sum of linear 
representations. The CG coefficients (we assume that all the CG coefficients are real numbers) are : 

(Al) 
(A2) 
(A3) 
(A4) 

where \m) = C^|a)|/3), \m) is the basis of a linear representation and |a), |/3) are the bases of two projective 
representations. 

In the following, we will illustrate the method to obtain the Hamiltonian of the tg phase as an example. 

The first step is obtaining the MPS. From Tab. II, the edge states of the to phase are described by the E3 projective 
representation. In an ideal MPS, every rung is represented by a direct product of two E3 projective representations, 
which can be reduced to four linear representations, £^3 (g) £^3 = Ag © Big © B2g © B^g. From Tab. Ill, Big, B2g,B-ig 
correspond to the bases respectively. The basis Ag (or |0, 0, )) is absent on every rung in the MPS 

state. Thus, the support space for the ideal MPS is the Hilbert subspace ®i{\l,z) © |l,y) © |l,a;))i, where i is the 
index of rung. From the CG coefficients (Al), we can write such an ideal MPS which is invariant (up to a phase) 
under the symmetry group: 

|^)= ^ Tr(A'"i...A™")|mi...mjv), 

{mi,...,mjv} 

with A'" = e*^'"i?C™. Here B is the CG coefficients of decomposing the product representations £3 ® E3 into a ID 
representation (here we choose B — C^^), and e*^"" can be absorbed into the spin bases. Now we have, 

A^^'y^^ay, A^^'^^^a,, (A5) 

The second step is constructing the parent Hamiltonian, which is a sum of projectors. Each projector is a projection 
onto the ground state subspace of two neighboring rungs. Assuming the orthonormal bases for the MPS state of two 
neighboring rungs i,i + 1 are ipi,ip2,'4'3,'4'Ai then the projector is Pi,i+i = ~(X]a=i \'^a) {'4'a\)i,i+i and the resultant 
parent Hamiltonian -ffoox — ^iPi.i+i is given as. 



5 2 2 

Y^(Sl,i + 82,1) • (Si^,;+i + S2,i+l) — 2^1,1 ■ ^2,1 — 3(^1>* ' ^2,i)(Sl,i+l • S2,i+l) 



■^(Sl^i • Si_i+i)(S2,i • S2,,;+l) + 3(^1,4 ' S2,j+l)(S2,i • Si^i+i) 



(A6) 



The final step is deforming the Hamiltonian. It can be shown that only the first two terms in (A6) are important. 
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To see this, we introduce the parameter d, 
5 



i 

2 



12 



(Si,i • S2,i+i + 82,1 • Si_i+i) 



^(Sl,^ • S2,i)(Si^,;+i • S2^i+l) + 3(^1,2 ' Sl,i+l)(S2,i • S2,i+l) + 3(^1,2 ' S2,i+l)(S2,,; • Si^i+i) 



(A7) 



Note that when d/J = 1, (A7) is the same as (A6). Now we study the ground state energy and entanglement 
spectrum (through time-evolving block decimation method) to see if there is a phase transition when d is varied. 




0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 




3E/3d 



0.4 0.6 

d/J 

(a) 




(b) 



FIG. 9: (Color online) (a) Energy VS. d/J. (b) The derivatives of the energy curve with respect to the parameter d and the 
information of entanglement spectrum. 



Fig. 9 shows that the energy is a smooth function of parameter d/J G [0,1]. Furthermore, the entanglement 
spectrum remains degenerate in d/J G [0, 1]. Thus we only need nearest-neighbor exchanges to realize the Haldane 
(to) phase, which leads to the Hamiltonian iJg in the main text. 

The active operators in Tab. II are obtained as the following. In the Hilbert space spanned by the two-fold 
degenerate edge states, only the Pauli matrices (cTa;, ctj,, CTz) can lift the degeneracy. But these operators are not 
physical quantities. We need to find physical spin operators which vary in the same way as these Pauli matrices under 
the symmetry group. In other words, we require that the active operators form the same linear representations as 
(fx, (Tj,, i7z), respectively. For example, in the projective representation, cFx varies as 

M(i?,)V,M(i?,) = a„ 
M(cr)V:,M(a) = (T^. 

On the other hand, from Tab. Ill, 

pt Z? C^: 

We find that and belong to the same linear representation B^g under the symmetry operation. This means 
that in the low energy limit (i.e. in the ground state subspace), these two operators have similar behavior. So we can 
identify as an active operator. Similarly, he operators and are active operators corresponding to ay and 
az respectively. 

Similar to (A6), we can construct the exactly solvable Hamiltonian iJ^cx of the phase. It is the same as (A6) 
except that every A • B term is replaced by AzBz — A^B^ — AyBy. This Hamiltonian can be simplified into the form 
of Hz [Eqn.(3) of the main text] without any phase transition (H^ and Hy are obtained similarly). 

The active operators in tz phase can be easily obtained: S^, S^, S^. Notice that S^, are not active operators, 
meaning that the edge states in the ground state will not respond to the uniform magnetic in x and y directions. To 
check this result, we perform a finite-size exact diagonalization of the solvable model Hzqx- As shown in Fig. 10, only 
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FIG. 10: The spectrum of the ground states and some excited states in tz phase, L is the length of the ladder, (a) B = Q, the 
ground states are 4-fold degenerate; (b) B \\ z, the degeneracy of ground states is split; (c) B A- z, the ground states remains 
degenerate. 



the magnetic field along z direction can split the ground state degeneracy. These properties are valid in the whole tz 
phase in thermodynamic limit. This verifies the conclusion that only is the active operator. 

This interesting result indicates that we can distinguish from <o by the response to magnetic fields. In phase, 
arbitrarily small magnetic field can split the degeneracy of the ground states, showing that the edge states carry free 
magnetic moments. According to Curie's law, the magnetic susceptibility will diverge at low temperature. But in tz 
phase, the edge states only carry magnetic moment in z direction, so the magnetic susceptibility within the XY plane 
does not diverge at low temperature, but it does diverge if the magnetic field is along z direction. These results are 
also verified numerically, see Fig. 11. The behavior of low-temperature magnetic susceptibility is measurable, which 




T/J T/J 

(a) (b) 
FIG. 11: The susceptibility in phase, (a) x diverges at T = if B |1 2; (b) x is finite at T = if S _L z. 
allows us to distinguish different SPT phases experimentally. 



Appendix B: Implementing the ladder Hamiltonian 

In this Appendix, we discuss in more detail the im- 
plementation scheme for the Hamiltonian with SPT 
phases. To keep our discussion general, we consider a 
two-dimensional (2D) coupled-harmonic-oscillator-array. 
Later, we will relate this general scheme to the specific 
example of a cavity-ladder as in the main text. 

Consider a 2D array, on each site of the array, two 
independent harmonic oscillators exist and couple with 



those on the neighboring sites via energy tunneling. We 
label these harmonic oscillators as L (longitudinal) and 
T (transverse), and assume that oscillators only cou- 
ple with those having the same label on the neighbor- 
ing sites along the direction specified by their labels. 
This can be achieved by requiring the frequency differ- 
ence between different types of oscillators (L and T) to 
be sufficiently large, and by setting up specific coupling 
schemes between neighboring sites. The Hamiltonian for 
this 2D coupled-harmonic-oscillator-array can be written 
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as {fi=l) 

j,k fi=L.T 

(Bl) 

where a^j, (o-Jk) the annihilation operator for the har- 
monic osciUator labeled L (T) on the jth site transver- 
sally and the fcth site longitudinally. The coupling 
strength along the longitudinal (transverse) direction is 
given by vl (vt)- 

The harmonic oscillators on each site interact with a 
two-level system {| tj,k}, \ ij.k)}, and the couphng rates 
are gj'/. — g^e'^^ ''* and gj^ — f;-^e'^3.'=*, respectively. 
This is illustrated in Fig. 12. The interaction Hamil- 
tonian is 

^2 - E E «k^^'''<A + H.c), (B2) 

j,k fj,=L,T 

where = | tj,k){ij,k |, and A-^ (A^) denotes the 
detuning of the corresponding harmonic oscillator mode 
(see Fig. 12). Finally, the two- level system on each site 
is coupled by a resonant dressing field, with the Hamil- 
tonian 

^3 = E E i^tkSt,+H.c.), (B3) 

j,fe ii=L,T 



where fJ^^j, = W^^e^'^^-'' is the Rabi frequency. 




I kk> ^ 

FIG. 12: Schematic for the coupling scheme of the two-level 
system on each site. 



Starting from the full Hamiltonian H = Hi+ H2 +-^3, 
we will eventually adiabatically eliminate the harmonic 
oscillator modes and derive an effective Hamiltonian for 
the dynamics of the coupled two-level systems through- 
out the array. 



Before doing so, let us first introduce the following transformations 



a' , = > exp -i^^ + a^„, (B4) 

|t;-,fe) = ^(e'^-Mt,,;c) + |;,,.)), (B5) 



where M and N are the total number of sites in the longitudinal and transverse direction, respectively. While Eq. 
(B4) diagonalizes Hi, (B5) and (B6) define the pseudo-spin basis || tjj.), \ 

With these, the Hamiltonians become 



=E E (B7) 

m,7i fj,=L,T 

^2-E E [^E-p"^^'^'''^^<ne^^^'(e-^"-^4^;rfc+e-""-^ (B8) 

j,k n=L.T vJWiV ^^^^ 

^a-E E ^^''S^k, (B9) 

j,k ti=L,T 



where ujh. 



basis states: S"^^ 



, and S'. 



2u-^ cos(^jp). The pseudo-spin operators are defined through the pseudo-spin 



3-k 
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We now go to the rotating frame via the transformation R = exp[— + ^3)^] 

_dR 



H" H'}j R-iR 



+ 1 

dt 



+ H.C.}. (BIO) 

Under the condition \ A'^ - A^| |A^| ± y^fl'^] ^ y^fl'^ > -^£^,a;^ „, the osciUator modes are virtually 

populated. We may adiabatically eliminate a^^^, and describe the dynamics of the system using the resultant effective 
Hamiltonian^^' 

^ ^ 2 AA" - c.^^,„ ^^■'■^'^^ ^ 4 AA- - c.#, „ + V2n^f ^ 4 AA" - c.f„ „ - Vin^f ^''^^''^^ 

J -ti fJi — -ij 5 TTL.TL ' ' 

AA' - ^fn.n + A/2r!''f ^ Ai^^Z^~7^' ^^3-k^r,k' 

I.!-! I I r ^ w'- c'+ 11 

4 ^ AA" - c^,'^,,^ - 7217^ AA" - a.,'^,,„ + ^n^^ ' ^-'^ 

(Bll) 

Adopting the formulae: = 0, = Z^n '^os(^)e^'^ = ^, and keeping only on-site and 

nearest-neighbor interactions, we can further simplify the effective Hamiltonian 

~ ^ ^ ^ 2[(AA')2- 2(17^)2] ^ J-'^-'' AA'[(AA^)2-2(f^'-f)2]^ ^■''^^ 2\/2[(AA')2 _ 2(17''f)2] 

V Z^l^ff ; 11 ^) (^ls,LY^hk'^3,k+l+ [(Ai)2 - 2(17''f)2j2 hk^ 3 ,k+l ^ ^ j 3,k+l) S 



3,k 



A^-LJ " 2^\9 J U--Lj " ^^T)2 '^J,*:'^i + l,fc + [(A^)2 - 2(f]''f)2]2 3,k^ 3 + ^,k ^ ^ 3 j + l.k) > ' 

(B12) 

While the first two terms on the first line in Eq. (B12) are constant and can be dropped, the third term is a Stark-shift 
in the pseudo-spin basis, and can be canceled via local optical elimination.'^^ This corresponds to applying a r.f. or 
Raman fields with appropriate magnitude and phase between the hyperfine states such that the effective Stark-shift 
is canceled. Then, under the condition 3(AA')2 _ 2[{V^Y , we have 

^off — [ i-) ' 2{Q^^y Z_/^^~ " ^3.k^j,k+l + '^j,k'^j,k+l + ^j.k'^jM+l)! 



jrk 

Eq. (B13) gives the most general form of the effective Hamiltonian using our setup. For the spin-ladder Hamiltonians 
we considered in the main text, we may take M = 2 so that j = 1,2 in the summations. Then, depending on the 
magnitudes and the relative phases of the coupling fields, we have either the Hamiltonian for the phase (io model) 

or the Hamiltonian for the tz phase {tz model). In particular, with — dj^k = ^J,k = Vj.k = for arbitrary 



\fj + l.k-Vj,k\ 
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{j, k}, the Hamiltonian reduces to the to model 



TT _ J \ ' \ ^ / C'Z QlZ I QlX QlX I Qly q'V \ 



= 1 k 



\\\^( C'Z QlZ I C'X QlX , Qly Q'y 



(B14) 



where the interaction rate J 



and A 



On the other hand, when 



■""i Wi.k ~ 'f2,k\ = a-nd iy9j^fc = ipj^k+i, the Hamihonian reduces to the model 









^j,k+i - ^j.fc 







TT _ t\ ^ \ ^( QlZ QlZ I QlX QlX , Qly Qly \ 

^Z — J 2_^Wj,k'^J,k + l + '^■J,k'^j,k+l + '^j,k'^J,k + l> 



3 = 1 k 



(B15) 



A straightforward example for the realization of the 
coupled-harmonic-oscillator-array is the coupled quan- 
tum electrodynamics (QED) cavity array, as shown in 
FIG. 5(a) in the main text. In such a system, atoms 
or solid spins interact with the quantized cavity fields, 
which couple to their neighboring ones across both the 
longitudinal and transverse directions via photon hop- 
ping.^* The two-level system in our general model can 
be replaced by a three-level structure, with two low-lying 
hyperfine states and an electronically excited state. Cor- 
respondingly, the coupling g^^, (.gj^,) in the general model 
is replaced by a Raman path in the longitudinal (trans- 
verse) direction, with an external laser field and a cavity 
mode each contributing a leg in the Raman coupling. 
Note that due to the large difference in the two-photon 



detuning of the Raman couplings, the two Raman paths 
are effectively independent. It is then straightforward 
to work out the correspondence: 5^ = ^ {-^ -I- ^), 

|A^| = - A^*!, where C is the Rabi frequency for 
the atom-cavity coupling, A^ is the detuning (c.f. Fig. 5 
in the main text). Importantly, one may realize Hamil- 
tonians of different SPT phases (to or t^) by adjust- 
ing the phases of the Rabi frequencies fi^^, and 
For typical experimental parameters'^": fi^ ^ lOOMHz, 
G'^ - lOOMHz, |Af| - IGHz {i = 1,2), - lOMHz, 
fi'f - lOOMHz, we have J ~ 0.15MHz, with the mag- 
nitude of A/ J widely tunable by adjusting the ratio be- 
tween VL^ and ri^. 
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